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We consider here chiral symmetry breaking in quantum chrome-dynamics arising 
from gluon condensates in vacuum. Through coherent states of gluons simulating a 
mean field type of approximation, we show that the off-shell gluon condensates of 
vacuum generate a mass-like contribution for the quarks, giving rise to chiral sym- 
metry breaking. We next note that spontaneous breaking of global chiral symmetry 
links the four component quark field operator to the pion wave function. This in 
turn yields many hadronic properties in the light quark sector in agreement with ex- 
periments, leading to the conclusion that low energy hadron properties are primarily 
driven by the vacuum structure of quantum chromo dynamics. 
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I. INTRODUCTION 



It is now accepted that quantum chromodynamics (QCD) is the correct theory for strong 
interaction physics of quarks and of hadrons. However, this being a nonabelian gauge theory 
of strong interactions, at present no reliable method other than lattice gauge theory |l ] 
is known for the solution of such problems ||. Although this is a beautiful and powerful 
method, the solutions finally need better computer capabilities compared to what is available 
now. It is thus desirable to consider alternative dynamical schemes which do not involve a 
discretisation of space-time. One such scheme with a limited objective has been the QCD 
sum rules |§ through duality principle. It optimises the use of a nonperturbative vacuum 
structure with perturbative QCD calculations |§, and reproduces nice experimental results. 
However, in a purely nonperturbative situation, the method is expected not to be applicable. 

We have proposed earlier an alternative scheme with a variational method which is 
nonperturbative @||, and uses off-shell gluons for the existence of a nontrivial vacuum 
structure |J. Calculations are done in Coulomb gauge. The inputs consist of the QCD 
Lagrangian, and QCD vacuum contains gluon condensates arising through a variational 
principle || . With a similar technique we have also considered chiral symmetry breaking in 
Nambu Jona Lasinio model, as well as for some quark antiquark phenomenological potentials 
used for hadron spectroscopy ]7||§. In the present paper we consider the same as arising 
from gluon condensates along with some low energy hadronic phenomena. 

The paper is organised as follows. In section 2 we recapitulate f| the nonperturbative 
gluon condensates in Coulomb gauge giving rise to a nontrivial vacuum structure and define 
notations. We then conjecture that through minimal coupling the quarks have a nonzero 
mass-like contribution, and obtain an approximate value for the same. In section 3 we 
derive the pion state as a quark antiquark pair through a constructional use of Goldstone 
theorem. In section 4 we proceed to use this relationship to obtain some low energy hadronic 
properties in the light quark sector. In section 5 we discuss the results. 

The method considered here is a non-perturbative one as we use only equal time quantum 
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algebra but is limited by the choice of ansatz functions in the calculations. The techniques 
have been applied earlier to solvable cases to examine its reliability || as well as to physically 
more relevant ground state structures in high energy physics |||9| and nuclear physics fllQjllH . 
Starting from the QCD Lagrangian we examine here the more complex problem of chiral 
symmetry breaking and its relationship to low energy hadronic properties. 



II. GLUON CONDENSATES AND CHIRAL SYMMETRY BREAKING 

In the present section we shall recapitulate the nontrivial vacuum structure in quantum 
chromodynamics through Bogoliubov transformation as considered earlier || to set the 
notations, and derive some results. 

The Lagrangian density including the quark fields is given as 

£ £ gauge ^matter ^int- (^T) 

Here, 

Cgauge = ~G a ^(d,W\ - d v W\ + g f ahc W\W°) + ^G\ U G^ . (2.2) 

where W a fl are the SU(3) colour gauge fields. Also the Lagrangian densities for the quark 
fields and interactions with gluons are given as 

^matter = V^7^ - 77l)V>, (2.3) 

and 

Ant = Q^in a W\i>. (2.4) 



In order to do the quantisation in Coulomb gauge, we write the electric field, G a oi in 
terms of the transverse and longitudinal parts as 

G a 0l = T G a 0t + dif a , (2.5) 



where the form of f a is to be determined. In the Coulomb gauge the subsidiary condition 
and the equal time algebra for the gauge fields are given as fll2"fl 



diW\ = 



(2.6) 



and 



W\{x, t), T G b oj (y, t)] = iS^Sa - -^)S(x - y). 



We take the field expansions for W a i and T G a oi at time t=0 as |13 

W a i{x) = (2tt)- 3/2 

and 



dk -> - 1 

(a a i(k) + a a i(-k) ) exp(ik ■ x) 

2u(k) 



T G a 0l {x) = (2vr)- 3 / 2 i J dk^\ ^-{-a\{k) + a^-kf) exp(z£ • x). 
From equation (2.7) these give the commutation relations for a a i and 



a a i(k),a b j (£j* 



S ab Aij(k)5(k - k ), 



where, ou(k) is arbitrary, and, 



A t3 (k) = fiy - hkj 



(2.7) 



(2.8a) 



(2.8b) 



(2.9) 



(2.10) 



13 k 2 

In Coulomb gauge, the expression for the Hamiltonian density, T 00 from equation (2.1) is 



given as 
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T 00 = : l - T G a jG a Ql + -W^-V 2 )^ + g f abc W\W c 3 d t W a , 



I 9 rabc 

+ ^{-i'fdi - m)%l) - gijft a W a iij : 



(2.11) 



where : : denotes the normal ordering with respect to the perturbative vacuum, say | 
vac >, defined through a a i{k) \ vac >= 0. The term \(dif a )(dif a ) automatically includes 
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interactions for both time-like and longitudinal gluons, here through the auxiliary field 
description, and T 00 is calculated after the elimination of the same 0. 

In Ref [3] the possibility of a nonperturbative vacuum having a lower energy was dis- 
cussed as from the above equation. Such a vacuum contained gluon condensates. We thus 
considered earlier a trial state, | vac > with such condensates over the perturbative vacuum, 
| vac > given as ||,|] 

| vac >= U | vac >, (2.12) 

where 

U = exp(fi t - B), (2.13) 

with 

fit f f(k)a a i (k)^a a i (-k) i dk. (2.14) 
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contains gluon pair creation operators with an ansatz function f{k) to be determined 
later. The operators, say, b a i(k), which annihilate | vac > are given as 

b a i(k) = Ua a l (k)U- 1 . (2.15) 

We can explicitly evaluate from the equations above the operators b a i(k) in terms of a a i(k). 
In fact, we then obtain the Bogoliubov transformations H] 



b a i(k) \ ( cosh/(A r ) -sinh/(fc)\ / a a i(k) 
b a i(-k)J V-sinh/(fc) cosh/(fc) / \a a t {-kf ' ' 



(2.16) 



where we have assumed that the function f(k) is even and real. Using equations (2.9) and 
(2.16), we obtain the same commutation relation for the operators b a i and b h given as 

^iik),^^') 1 =5 ab A lj {k)5(k-k'), (2.17) 

which merely reflects that the Bogoliubov transformation (2.16) is a canonical transforma- 
tion. In order to consider the stability of | vac > with respect to the transformation of 



equation (2.12), the expectation value of T 00 with respect to | vac' > was evaluated and was 
then minimised over f(k). 

It was noticed || that under certain conditions the perturbative vacuum becomes un- 
stable, where it was also seen that the general solution of f(k) through minimisation is 
impossible to obtain analytically. We had taken an ansatz function f(k) of equation (2.14) 
for | vac' > description as given by || 

sinh(f(K)) = Ac- Bk2/2 (2.18) 

and then had determined the parameter A through a minimisation of vacuum energy density. 
The dimensional parameter B was determined through the evaluation of the SVZ parameter 
||. We may note from equations (2.16) and (2.18) that 

< vac\a a i(k) ] a a i(k')\vac' >= 16. A 2 exp{-BK 2 ) 6 (k - k'), (2.19) 

so that the ansatz of equation (2.18) merely implies a Gaussian distribution for the per- 
turbative gluons in \vac' >. The minimisation was done for pure gluon fields in Ref. || 
with a s = g 2 /(4n) of the Lagrangian as the input coupling strength. For example, we 
then have for input values of a s = 0.5, 0.6, 0.8, 1.0, the identification for the output 
parameters as A min = .178, .306, .557, .658, B^ = .40, .48, .62, .69 in fermis, and 
m = 232, 272, 325, 345 in MeV respectively, where uj(k) = \J k 2 + m 2 . "m" here is a mass- 
like parameter for uj(k) of equations (2.8) which was determined through a self-consistent 
calculation. It is not associated with energy momentum four vector of a "free" gluon, which 
anyhow does not exist. We may note that the results of Ref. || are very similar to some of 
the earlier results and that the method of solution for the auxiliary equation avoided 
Gribov ambiguity [ |I~6|j . 

Here we shall be using the final output, the stable vacuum | vac' > as the physical 
vacuum. We shall assume that the vacuum structure in QCD in presence of quarks is 
mainly driven by gluon condensates, and study some consequences of the same. 

We now proceed to consider the quark fields in the presence of the "dressing" of vacuum 
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with off-mass shell gluon quanta as considered above. We note that here with minimal 
coupling 

P -> P+ < vad\gH a t h W a i W h i \vad > . (2.20) 

With the mean field approximation for the gluon fields as above, we then see that the mass 
of the quarks mq is now given as 

m 2 Q ~m 2 + g\t a f u (0), (2.21) 

where we have substituted || 

< vac' | W\(x)W h 3 (y) \ vac' > = — - J dkd k ^ A^k) 



In the above 



(2tt) 3 J oj{k) 
5 ab Ux-y). (2.22) 



F + ik) = ^^ + s^f(k)). (2.23) 

With t a t a = 4/3, we then obtain that 

m| = m 2 + ^pM0) (2-24) 

where summation over the index % is understood. We thus note that even when m = there 
is a finite mass of the quark, which implies breaking of chiral symmetry. We then obtain 
through numerical evaluation || of fa(0) that m 2 Q = .0358 GeV 2 for a s = 0.5, so that we 
then have 

m Q = 189MeV. (2.25) 

We may however note that the above result is deceptive. It is incomplete since we 
have made the mean field approximation with minimal coupling to generate a mass-like 
contribution, and do not know how it operates regarding details of dynamics. We use this 
only to illustrate that gluon condensates break chiral symmetry. We shall include this effect 
in the expansion of the quark field operators. 



We then write the quark fields as 

1 



ip(x) 



(2tt) 3 / 2 



U(k)qi(k) + V(—k)qj(—k) 



(2.26) 



where 17,18 



U{k) 



cos 



o ■ ksin 



(2.27a) 



and 



V(-k) 



—a ■ ksin^Y- 



cos 



X(k) 



(2.27b) 



The above form is so taken that it satisfies the equal time algebra [17, IB]. In the above, for 
free quark fields of mass m<g, we have [ |T7| ] 



cos- 



(2.28) 



2 v 2j9° ' 1 2 v 2p° 

with p° = (nig + k 2 ) 1 . We may further note that when chiral symmetry is unbroken x(k) = 
7r/2. However, since the quarks are not free, and, chiral symmetry is broken, x{k) shall be 
different from the above. 

We would like to associate the quark field operators as in equation (2.26) explicitly with 
the presence of quark condensates generating chiral symmetry breaking. Parallel to the 
earlier operators for gluon condensates, let us consider the quark antiquark pair creation 
operator, with h(k) trial function, 



B ] Q = J h(k) qi (kya ■ kq!(-k)dk. 



(2.29) 



Let us now consider the unitary transformation 



U, 



Q= e 



BL—Br 



(2.30) 



which will operate on chiral vacuum. Also, let us write the corresponding quark field operator 
of equation (2.26) in momentum space as 



Mk) 



^ cos^-qj(k) — a ■ ksin^ qi (—k) ^ 



a ■ ksin^ qi {k) + cos^-qi{—k) 



(2.31) 



When chiral symmetry is there, xo — 7r/2. We may now note that the unitary transformation 
Uq is equivalent to the Bogoliubov transformation given as 



( 



qi(k) 
qi(-k) 



cos(h(k)) a • ksin(h(k)) 
—a ■ ksin{h{k)) cos(h(k)) 



^ ( qi(k) ^ 



qi(-k) 



(2.32) 



Now, with ip(k) = Uq^o^Uq, we generate equation (2.26) when we identify that 

2 2 v ' 



(2.33) 



We thus note that the form of equation (2.26) with (2.27) corresponding to chiral symmetry 
breaking can be interpreted as a destabilisation of chiral symmetric vacuum through equation 
(2.30) when xo — tt/2. Hence taking the quark field operators with chiral symmetry breaking 
or taking a condensate over chiral vacuum become equivalent. The quark field operators 
after chiral symmetry breaking get related to the vacuum structure. 
We may also note that 
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< ip(x)ip(y) >=< vac\U^ip (x)il;o(y}UQ\vac > 



cosx(k)e lk ^-^>dk, (2.34) 



where, the factor 12 arises from colour, flavour and spin degress of freedom. We may 
substitute cosx{k) = sin(2h(k)), describing the above function explicitly in terms of the 
quark antiquark correlations. We also have 



< ip(x)ip(x) >= — /i 3 



(2.35) 



with fj, determined in terms of the condensate function or cosx{k). The condensate function 
h(k) vanishes when cosx{k) = 0, as should be the case for chiral vacuum. 

The form of x{k) will, generally speaking, depend on interactions [0,0], and, can be 
obtained through an extremisation |pl,|T9"|l. Also, the function x{k) will get related to the 
pion wave function using Goldstone theorem as shown below. This shall include vacuum 
structure as a post facto information linked to phenomenology without an extremisation of 
energy density containing highly nonlinear expressions ||. 
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III. CHIRAL SYMMETRY BREAKING AND PION WAVE FUNCTION 



We shall here first consider only two quark flavours and, the usual three colours. When 
chiral symmetry remains good, we note that 

Q 5 i \vac>=0 (3.1) 

where Q§ is the chiral charge operator given as 

Q 5 l = J ^ T -^{x)dx. (3.2) 

In this case quarks are massless and we have 

cos— = stn— = -j=. (3.3) 

For symmetry broken case however [|7|,|8| 

Qj | vac >^ 0. (3.4) 

This will correspond to the pion state. To show this we first note that 

Q 5 \H]=0. (3.5) 

Clearly, for the symmetry breaking phase, \vac' > is an approximate eigenstate of the total 
Hamiltonian H with EV as the approximate eigenvalue, V being the total volume and E the 
energy density ||. With H e ff = H — EV, we then obtain from equation (3.5) that 

H eff Q 5 i \vac' >= (3.6) 

i.e. the state Q§ \ vac' > with zero momentum has also zero energy, which corresponds to 
the massless pion |7|,§]. Explicitly using equations (2.26) and (2.27), and with qi now as two 
component isospin doublet corresponding to (u,d) quarks above, we then obtain the pion 
states of zero momentum as 

| 7^(0) >= N n -^= f q I {kyr i q I (-k)cos{ X {k))dk | vac >, (3.7) 



\/6 
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where, N n is a normalisation constant. The wave function for pion u n (k) thus is given as 
proportional to cosx(k). The colour, isospin and spin indices of q^ and q for quarks and 
antiquarks have been suppressed. Further, with 

< n\0) | 7T j (p) >= 5 ij 5{p), (3.8) 



the normalisation constant N n is given through [17 



N v 2 cos 2 {x(k))dk = l. (3.9) 



In the notations of Ref. [17j thus the pion wave function u n (k) is given as 

u n (k) = N n cos( X (k)). (3.10) 

Clearly the state as in equation (3.7) as the Goldstone mode will be accurate to the 
extent we determine the vacuum structure sufficiently accurately through variational or any 
other method so that \vac' > is an eigenstate of the total Hamiltonian. The above results 
yield pion wave function from the vacuum structure for any example of chiral symmetry 
breaking, and is a new feature of looking at phase transition through vacuum realignment 



We had earlier considered low energy hadronic properties [[17]] with an assumed form of 
interacting quark field operators. This is equivalent to a choice of cos (%(&)). In addition, 
we had taken Gaussian wave functions for baryons and mesons, and then discussed the 
phenomenology. We shall here examine the consistency of the present picture with low energy 
hadron phenomenology where pion wave function and quark field operators get related as 
in equations (2.27) and (3.10). 

The notations of Ref. ||17|| , which we shall use now, correspond to 



f q (k)=cos^-; \k\g q {k) = sin^-. (3.11) 

We may note that as in equation (2.26) with (2.27), we are retaining the fully relativistic 
four- component quark field operator. The phenomenological assumption here will consist of 
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explicitly taking a specific expression for cosx(k) . We note that cosx(k) as in equation (2.34) 
is the correlation function for quark pairs in vacuum. Hence, parallel to the choice of gluon 
correlation functions in equation (2.19) as a Gaussian ||, we shall take here the ansatz 

cos X (k) = exp (^-^fk^J . (3.12) 

It also implies through equation (3.10) that the pion wave function is a Gaussian. From 
equation (3.9), here we have 

R 3/2 

N^ = : ^ JI = 0A2AxR 3 J 2 . (3.13) 

From equations (2.34) and (2.35) we also obtain that 

(12)V3 



R- 1 = 0.913 x R- 1 . (3.14) 



V2n * 

We would now like to see what may be the nature of dispersion curves for quarks in 
vacuum. We note that for free quark fields sinx(k) = n/e{k). We may use this to obtain 
that 

e(k) = K = , (3.15) 

Sjl - cos 2 x(k) 

giving the dispersion curve for quarks in vacuum as a medium. Thus for equation (3.12), 

t(k) = K = . (3.16a) 
V 1 — e R * k2 

Hence, for small k 2 , 

e(£) = i- + ^ 2 . (3.16b) 

— # 

We thus see that e(k) does not have the form of a free particle. We note that for large ac, 
e(k) ~ k, which follows for any choice of cosx(k) from equation (3.9). 

We shall need later expansion of cosx{k + p) for small \p\ for charge radius and magnetic 
moment. For this purpose we use the notations 

diCosx(k) = —ki bi(k); didjCosx(k) = —5ij bi{k) + kikjb 2 (k). (3-17) 
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In the above, we have used rotational invariance. We also note that 

k 2 h(k) 2 = (ycos X (k))\ (3.18) 

and, define c 2 {k) such that 

h(k) - l -k 2 b 2 (k) = c 2 (k) = - l -^ 2 cos X (k). (3.19) 

The above equations shall be useful if we wish to change the ansatz function for cosx{k). 
Clearly, for equation (3.12), we have 

h = Rlcos X (k); b 2 = -R*cos X (k). (3.20) 



IV. LOW ENERGY HADRONIC PROPERTIES 



In the context of chiral symmetry breaking, the model of Ref. |L7| now has got modified 
since the pion wave function and the quark field operators in the light quark sector get 
related. This in many ways enriches the earlier model [[17[], as it decreases the number of 
independent quantities. The previous model had also been extended to incoherent processes 



involving structure function as well as hadronisation p0|j21[] . It has also been applied to the 
problem of strong CP violation recently p2l . We here consider different low energy hadronic 
properties parallel to JT7| using equations (3.11). 



A. Pion decay constant 

The pion decay constant c n as calculated earlier in terms of the wave function with 
relativistic correction to Van Royen-Weisskopf relation |23] is given by [ |H]j 



;i + 2nVK(o) 



y/E 



(4.1) 



where the value of the wave function is taken at the space origin, and, g q is a differentiation 
operator. In fact taking the Fourier transform the left hand side of equation (4.1) becomes 
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1 -2g q (k) 2 k 2 )u n (k)dk. 



(27r) 3 /2 

From equation (3.11) however 1 — 2g q (k) 2 k 2 = cos(x{k)). Hence from equation (3.9) the 



normalisation constant N w is given through [17], with c n = 92 MeV, 



(2tt) 3 / 2 N n 
so that in equation (3.10) 



c„{m T ) 1/2 



0.0140GeV 3/2 , 



(4.2) 



N v = 4.55 GeV- 3/2 . 



(4.3) 



Hence, from equations (3.13) and (3.14) we may note that, for the choice of equation (3.12), 



Rl = 7T x iV^ 3 = 23.69 GeV" 2 ; y, = 187 MeV. 



(4.4) 



We may note that in equation (3.16b) if we take rriQ = then mq for the above 

equation is around 205 MeV, which may be compared with the earlier value of 189 MeV 
from the vacuum structure in QCD due to gluon condensates as in equation (2.25). 

Thus the known value of the pion decay constant puts a severe constraint both on the 
pion wave function as well as on the form of the quark field operators. 

B. Pion charge radius 



We note that the pion form factor is given by JT7 

1 m 



<n+(-p)\J°(0)\n+(p)>- 



(2vr) 3 p° 



TGlit) 



where through direct evaluation [17 we obtain that 



In the above [I7| , 



t = 



«i = h ttP- 



(4.5) 



(4.6) 



(4.7) 
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We shall here retain only terms upto \p\ 2 and use 

G E (t) = 1 + . 

For the expansion in powers of t, let us substitute in equation (4.6) 

- 1 - 1 

k i = k ~ 2^ h = k + -p 

corresponding to equation (4.7). Using equations (3.17) and (3.19), we get 
cosxik'^cosxiki) = cos 2 x(k) — ( cosx(k)c2(k) + -k 2 b\{k) 2 



4 V /vw v ' 3 
and, 

cosx(^i) + cosx(A;i) = 2cosx(^) -^-c 2 {k). 

We then have in the lowest order in 



u^KVuM = n 2 



cos 2 X (k) - i^- (co SX (£)c 2 (£) + -fc 2 fei(fc) 2 



Also, by equations (3.11) and (4.9), 

fQ(ki)fQ(ki) = — c 2 {k) + 



2 16 ^ 3(1 + cos X (k)) 

— * — * 

with k[ = \k[\ and K\ = \ki\, 

<9 Q {k M k) = - f U ( £) + Ph ^ )2 



2 16 \ 3(1 -cosx{k)) 

and, 

fci . jfei = 1 - 

6k 2 

We then easily obtain from equations (4.6), (4.7) and (4.8) that 



p2 p2 I p2 



where, using equations (4.12), (4.13), (4.14) and (4.15) and simplifying, 
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N 2 

R i = ^T I (Vcos X (k)) 2 dk 



is the contribution coming from the wave function alone, and, 



AT 2 

R\ = I cos 2 X (h ) 



k%(k) 2 2(1 - cos X (k)) 



1 — cos 2 x(k) 



k 2 



dk 



is the balance of the contribution. 

For the choice of equation (3.12) we then have 



(4.17) 



(4.18) 



R lh 



(4.19) 



where I = .0159 as numerically evaluated. This yields that R 2 h = 9.26 GeV 2 , or, R ch = .605 
fms, which may be compared with the experimental value of R c h = 0.66 fms or, R 2 h = 
11.22 GeV" 2 M. 



C |gA/gy| 

We next consider some properties associated with proton and neutron. Let us approxi- 



mate the proton state with a harmonic oscillator wave function, so that we take [17| 



/3i? 4 \ 3//4 / R 2 
u p {h, k 2 , h) = l^-fj exp l--* J2(h - k 3 ) 2 ] . (4.20) 



With the relativistic corrections, we then obtain that [17 



\9a/9v\ = I ■ (^) 3/2 / exp (~^ 2 ) Q + \cosxik)) dk. (4.21) 

We have obtained here the above expression through a direct evaluation of 

<Pi/ 2 (0)|J 5 >i/ 2 (0) > 

retaining the four component quark fields in the current through equation (2.27), and have 
substituted from equation (3.11) that f 2 — \gqk 2 = |(1 + 2cos (%(£;))). Taking f q =l (and 
hence g q = 0) or, cosx{k) = 1, will give the result for nonrelativistic quark model. Here the 
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full relativistic structure of the quark field operator is being retained. \qa/ gv\ = 1-257 puts 
a constraint on the quark field operators and the approximate size of the nucleons as below. 

In fact, we can use the above to obtain the value of R 2 for the choice of equation (3.12). 
We note that we then have 

3/2" 



1.257 = - 
9 



1 + 2 



3a 



,3a + 1, 

where we have taken R 2 p = aR\. This yields that R 2 p = .5707 x Rl = 13.52 GeV~ 2 . 



(4.22) 



D. Proton charge radius 

In order to make another use of the size of the proton, let us calculate its charge radius 
with relativistic corrections. The charge radius is given through the same identification as 
in equation (4.8) with (17| 

1 m r 



< P r(-p)\r(0)\Ps(p)>=Sr. 



(2tt) 3 p 

In the above, a direct evaluation yields that, upto order \p\ 2 



fGUt). 



G p E (t) 



'3# 
~2T 



2\ 3/2 



exp 



o p2 



4 



x 



f q (K)f q {h) + K ■ hg.ik'Mh) 



dk. 



(4.23) 



(4.24) 



Here, in contrast to equation (4.9), we have ]T7|] 



- 2 - 2 

k i = k - -p; h = k+ -p. 



We then easily see as before that the proton charge radius is given as 



(4.25) 



p2 _ p2 , r?2 



(4.26) 



where R\ = R 2 is the nonrelativistic contribution from the proton wave function as may be 
seen with a contribution parallel to that of equation (4.17), and, 
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Rn 



1 {3R; 



2\ 3/2 



X 



9 V 2tt J 
k 2 b x {kf 



exp 



1 — cos 2 x(k) 
is the "relativistic" contribution. 



+ 



2(1 - cos X {k)) 
k 2 



dk 



(4.27) 



We have determined R 2 in the previous subsection. Substituting the same and using the 
form of cosx{k) as earlier, we obtain through a numerical evaluation that R\ = 3.08 GeV~ 2 , 
so that we obtain R 2 ch = 16.6 GeV -2 . This gives that 



Rch = -81 fms., 



(4.28) 



which is similar to the experimental value of the same. 



E. Proton magnetic moment 



We identify the magnetic moment fi p of the proton through the equation [I7| 



a x (-2£)} rs n P =< Pri-pjlA^lPsip) >, (4.29) 



(27T, 

where on the right hand side above all terms involving \p\ 2 will be neglected. Such an 
evaluation in fact through the notations of Ref. []T7|| yields that 



'3R 2 \ 3/2 r ( 3i? 2 



a x {-m n lh = h7 / expl-^k 2 ) dh^Xk'M^fy. (4.30) 



In the above, the momenta k\ and k[ are given by equatin (4.25) as for the charge radius, 
and, u L ( p \k) = S(L{p))u(k) includes spin rotations of the quarks [|17|]. On carrying out the 
simplification of the above expression with a straightforward but complicated algebra, we 
obtain that 

Atp = n P i + fip2, (4.31) 

where 
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and, 

In the above we have used equation (4.13) and (4.14) as earlier. We may note that [i v \ above 
came from Lorentz boosting and in the nonrelativistic limit equals one nuclear magneton. 
fi P 2 in the nonrelativisic limit equals (2/3) of quark magnetic moment as in the earlier 
calculation [17]. We may note that 

With the evaluation of the integral in equation (4.33) using equation (3.12), we obtain from 
equation (4.31) that 

H P = — x 2.88. (4.35) 
2m 



and, from symmetry or through direct evaluation [17 



H n = x -1.92. (4.36) 
2m 

These results are higher than expected, but still in reasonable agreement with experiments. 



V. DISCUSSIONS 



We should note that the present approach to chiral symmetry breaking through vacuum 



destabilisation []7||| is different from the earlier attempts ||25|| , which mainly use Schwinger 
Dyson equations to get the gap equation, and current algebra to obtain some additional 
results. We have more conclusions because of the explicit construction of vacuum state, 
here even relating the low energy properties of hadrons to a distribution function for quarks 
and antiquarks in vacuum! 
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Our results here are only linked to the vacuum structure. There will be additional 
dressing of the hadrons by gluons P,|26|| which will affect the wave function of the hadrons. 
Chiral symmetry breaking also is only approximate, which will change the results for pion 
||. Also, old SU(6) symmetry in spin and flavour space has been used for the construction 
of proton and neutron states, where we know that this symmetry is not sufficiently good. 
This shall have additional effects not included here ||17| , |2~0 . 



Here cosx{k) through equations (2.33) and (2.34) is really governed by the vacuum 
structure of QCD for chiral symmetry breaking. The surprising feature is that this structure 
yields many hadronic properties in the light quark sector in agreement with experiments, 
leading to the conclusion that the low energy hadronic properties are primarily driven by 
the vacuum structure of quantum chromodynamics! 

The results suggest that we should try to get cosx(k) through a variational procedure 
by minimising energy density parallel to ||. Such a programme would relate R n to the 
vacuum structure of QCD where the value of R n will depend on a s ||. Implementation 
of the same will enable us to recognise the coupling constant a s for which it agrees with 
present phenomenological value, which will be an additional conclusion from the present 
calculations. 
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